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The observation of the v e energy spectrum from a supernova burst can provide constraints on 
neutrino oscillations. We derive formulas for adiabatic oscillations of supernova antineutrinos for 
a variety of 3 and 4-neutrino mixing schemes and mass hierarchies which are consistent with the 
LSND evidence for oscillations. Finally, we explore the constraints on these models and 

LSND given by the supernova SN1987A P e ’s observed by the Kamiokande-2 and IMB-3 detectors. 
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I. INTRODUCTION 


In recent years, the treatment of neutrino transport in 
the environment of a core-collapse supernova (SN) explo¬ 
sion has improved to the point of making realistic predic¬ 
tions on the observables for neutrinos reaching the Earth 
§110- Of particular interest for this paper are the 
average energies at the neutrinospheres, i. e. the surfaces 
of last scattering for the neutrinos, estimated to be 10—13 
MeV for v e , 14—17 MeV for v e , 23—27 MeV for Ucr 


The differences in temperatures between the various 
neutrino flavors can be qualitatively understood. Heavy- 
lepton neutrinos can interact only via neutral-current 
(NC) processes, the main contribution to their transport 
opacity coming from neutrino-nucleon scattering, which 
dominates over neutrino-electron scattering. In addition 
to this same NC contribution, the transport opacity for 
fe's and F e ’s depends also on the charged-current (CC) 
absorptions v e + n —> p + e~ and v e + p —> n + e + , re¬ 
spectively. Therefore, the v e - and T^-spheres are located 
at larger radii with respect to the other neutrinospheres, 
that is at lower densities and lower temperatures. More¬ 
over, in a neutron-rich environment, v e + n —> p + e~ 
dominates over v e + p —> n + e + : the emergent u e ’s orig¬ 
inate from layers further outside the center of the star 
compared to v e 's, therefore at lower temperatures. The 
total energy released in a SN explosion is approximately 
equipartitioned between the different neutrino and an¬ 
tineutrino flavors ||. 

The above predictions can be confronted with the ob¬ 
servation of the supernova u e energy spectrum detected 
on Earth. Neutrino oscillations are expected to modify 
the spectrum since (E Ce ) < The energy depen¬ 

dence of the neutrino cross-section in the detector mate¬ 
rial, approximately ap eP oc (Ep e — 1.29MeV) 2 |J, helps 
in making the v e energy spectrum distortion a sensitive 
experimental probe to neutrino oscillations. This is be¬ 
cause higher energy neutrinos interact significantly more 
than lower energy ones. 


We show that the extent of the spectrum modification 
depends crucially on the specifics of the neutrino mixing 
scheme and on the neutrino mass hierarchy under con¬ 
sideration, and we derive the relevant formulas assum¬ 
ing an adiabatic propagation for the antineutrinos in the 
supernova environment. Antineutrinos propagate adia- 
batically if the varying matter density they encounter 
changes slowly enough so that transitions between lo¬ 
cal (instantaneous) Hamiltonian eigenstates can be ne¬ 
glected throughout the entire antineutrino propagation. 

So far, neutrinos from one supernova have been de¬ 
tected and their energy measured: SN1987A was ob¬ 
served by the Kamiokande-2 and IMB-3 detectors. The 
overall 20 events seen by those two detectors have all 
been interpreted as v e interactions |J. We examine the 
constraint of such observations on the LSND allowed re¬ 
gion of > v e oscillations {tJ , for various neutrino mass 
and mixing models. If the LSND evidence is confirmed by 
the MiniBooNE experiment || , several models can be ex¬ 
cluded or constrained on the basis of the observations of 
the supernova SN1987A and possibly future supernovae. 


II. ADIABATIC OSCILLATIONS AND 
NEUTRINO MIXING SCHEMES 


A. energy spectrum and the permutation factor 


In the presence of neutrino oscillations, the v e flux 
reaching the Earth, F „ e , can be different from the pri¬ 
mary flux at the neutrinosphere, . We will assume 
that, at production, the energy of active antineutri¬ 
nos is equally divided into the three active flavors, i.e. 
that J 0 °° dEjj a Ej; a F® a has the same numerical value for 
a = e, /i, t. Moreover, we will also consider neutrino 
mixing models where the three active neutrino species 
are augmented by a fourth sterile neutrino with no stan¬ 
dard weak couplings: in those cases, we will assume that 
the sterile component is negligible at production. 

The neutrino flux reaching the Earth is: 
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Fv e ~ {Pfi—^e + Pr —+ Pe— ►e-^p e 

oc(pi^+(l-p)i£) 
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where we have defined the permutation factor p as: 


P p —*e T Pr —>e H“ Pe —>e 

and p^^T^e-te are the probabilities for a i+, u T , v e respec¬ 
tively at the neutrinosphere to oscillate into a v e . In 
Eqs.|l],[|, we have assumed that p is energy-independent 
(as will be justified later), and that {Ep } = {Ep r ). In 
Eq. |l], we neglect the (energy-independent) proportion¬ 
ality factor since we will not deal with event rates, but 
only with neutrino energy distributions. 

B. Neutrino propagation in the adiabatic 
approximation 

In vacuum, the Hamiltonian that governs neutrino 
propagation is diagonal in the mass eigenstate basis | vf): 

( H 0 )ij = {ui\H 0 \i/j} = Eidij ( 3 ) 

If the neutrinos all have the same relativistic momen¬ 
tum p , their energies Ei differ only by a term proportional 
to their squared-mass differences, since Ei ~ p + mf/ 2 p. 
If U is the unitary mixing matrix that relates the fla¬ 
vor eigenstates \v a ) to the mass eigenstates via \v a ) = 
U a i\i>i), the elements of the vacuum Hamiltonian in the 
flavor basis are given by |J: 

2 

777 t 

(Ho)af> = KiUpi-^ ( 4 ) 

where we have neglected the contribution pS a0 in (Ho) a0 , 
which is irrelevant for neutrino oscillations. 

In matter, D e ’s undergo coherent CC forward¬ 
scattering from electrons, and all active flavor antineutri¬ 
nos coherent NC forward-scattering from electrons, pro¬ 
tons, and neutrons in the medium. These processes give 
rise to an interaction potential V = Vw+Vz, which is di¬ 
agonal in the flavor basis and proportional to the matter 
density p: 

(V) a0 = A a ^5 a0 ( 5 ) 

tun 

where A a is a proportionality constant, in general differ¬ 
ent for a = e, p, r, or s, Gf the Fermi constant, and to at 
the nucleon mass. The relevant Hamiltonian for neutrino 
propagation in matter is therefore H = Ho + V. 

At the neutrinosphere, the density p is so high (~ 
10 12 g/cm 3 0]) that the interaction potential dominates 
over the vacuum Hamiltonian, so that the propagation 
eigenstates coincide with the flavor eigenstates. As the 
propagation eigenstates free-stream outwards, toward re¬ 
gions of lower density, their flavor composition changes, 
ultimately reaching the flavor composition of the mass 
eigenstates in the vacuum. Given that the neutrinos es¬ 
cape the SN as mass eigenstates, no further flavor oscil¬ 
lations occur on their path to the Earth. 


Model 

Hierarchy 

Propagation 

Normal (1+1+1) 

m 3 > m 2 > mi 

—► v 3 

vp —► v 2 

v a —► id 

Normal (1+1) 

m2 > mi 

—> V2 

v e —> O! 

LSND-inverted (1+1) 

mi » m 2 

—> Vl 

V e -*• v 2 

Normal (2+1) 

m 3 > m2 » mi 

tV -*• tb 
v T —► v 2 

—► Vl 

LSND-inverted (2+1) 

mi > m 3 > m 2 

—*■ t'l 

Vt —► V 3 

Ve —► V2 

Normal (2+2) 

m 3 > m 2 > mi > mo 

Vp -> v 3 

V T —> V2 

V a —► V\ 

Ve —> Vq 

LSND-inverted (2+2) 

mi > m 0 > m 3 > m 2 

Vfs —> Vi 

V T —► Vq 

V s —>■ V 3 

Ve —► V2 

Normal (3+1) 

rri4 +> m 3 > m2 > mi 

Vfs —> V A 

Vt —*■ V 3 

V B V 2 

Ve —> Vl 

LSND-inverted (3+1) 

m 3 > m 2 > mi m 4 

Vp —► V3 

V T —► V2 

Vs —► Vl 

V e —► V4 

TABLE I: Adiabatic neutrino propagation in 
for the neutrino mixing models considered. 

the SN ejecta 


More specifically, making use of the adiabatic approxi¬ 
mation and of the fact that no energy-level crossing is 
permitted, the flavor eigenstate at the neutrinosphere 
with the maximum interaction potential reaches Earth 
as the mass eigenstate with the biggest neutrino mass. 
In general, the energy level order is maintained through¬ 
out the neutrino propagation in the SN ejecta. 

This is illustrated in Tab. | for three neutrinos in the 
row labelled “Normal (1 + 1 + 1)”, where we have taken 
A 1 > Ap > A a and m 3 > m 2 > m\. For example, the 
probability for a v a to emerge from the SN environment 
as a v 0 is given by: 


p a ^0 = \(9p\U evol \9 a )\ 2 = \(Upm\U ev ° l \9 a }\ 2 = 

\UpAi\ 2 = \U 0 i\ 2 ( 6 ) 
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where U evo1 is the adiabatic evolution operator. In Eq]t], 
we have used Tab. | to get: 

(Di\U evol \u a ) = S iA (7) 

This result can be immediately generalized to any 
number of antineutrino generations. Also, as long as the 
adiabatic approximation is satisfied, the formula does not 
depend on the specific dynamics for the neutrino prop¬ 
agation, for example on the number and position in the 
SN environment of MSW-resonances. We will comment 
more on the validity of the adiabatic approximation in 
the next section. 

In this paper, we consider three or four flavor com¬ 
ponents, including a sterile one. At tree-level, the pro¬ 
portionality factors A a in the interaction potential for 
neutral matter are 1-0: 


I (1 - 3Y' e )/v / 2, for & e 

^ 4 =< (l-Y e )/y/2, forPy,P T ( 8 ) 

[ 0, for D s 

where Y e is the electron fraction per nucleon. Follow¬ 
ing the assumptions of 0 0, we use Y e ~ (1 + 
(F Pe )/(£ , [/e )) _1 > 1/3 at the neutrinosphere. Consid¬ 


ering also one-loop electroweak radiative corrections, a 
difference in the Py and v T interaction potentials of mag¬ 
nitude (Ay — Aj) IAn ~ 10~ 4 appears due to the dif¬ 
ference in the charged lepton masses 0 0- At the 
neutrinosphere, this second-order effect in the interaction 
potential dominates over the vacuum Hamiltonian terms 
(as long as \m 2 — ra 2 \ < 10 eV 2 for all i,j), and removes 
the Py — P T degeneracy. Therefore, for the antineutrino 
channel considered here, we take: 

Ay > A T > A s > A e (9) 

For the neutrino channel, one should substitute A —> — A 
in Eq.|[ and the order in Eq.|] would be inverted. 

Therefore, given a specific neutrino mass and mixing 
model, the permutation factor can be easily evaluated 
in the adiabatic approximation, and its numerical value 
does not depend on the neutrino energy. We will com¬ 
ment on possible energy-dependent Earth matter effects 
in the next section. In practice, one proceeds backwards: 
given a certain measured value of p, it is possible to con¬ 
strain possible models for neutrino oscillations. This ap¬ 
proach is used for example in J0] to constrain models 
explaining the solar and atmospheric neutrino data; in 
this paper, we focus on 3 and 4-neutrino models explain¬ 
ing the LSND data. 


C. Possible mixing schemes 


parameters and p for the eight possible mass and mix¬ 
ing schemes considered below are given in Tab. [n| The 
mass hierarchy and the adiabatic propagation of the neu¬ 
trino eigenstates for these mixing schemes are depicted 
in Tab. |. 

The simplest possible mixing scheme is a (1 + 1) model 
explaining only Py —> v e LSND oscillations in vacuum, 
and not the atmospheric or solar oscillations: 



cos ■& sin -d 
— sin •& cos •& 



( 10 ) 


where the mixing angle •& can assume any value in the 
range 0 < § < 7r/4. 

We consider a (2 + 1) model motivated, for example, by 
CPT-violating scenarios (see, e.g. [ 14 ), in which atmo¬ 


1 N 



= 

\Vr ) 

\ 



spheric and LSND oscillations in the antineutrino chan¬ 
nel are obtained via the mixing []l5f : 


( 11 ) 


The matrix in EqJTl] is chosen to ensure large Vy —> v T 
mixing for atmospheric neutrinos (sin 2 2'd atm = 3/4), 
while the LSND Py —► P e mixing is fixed by the parame¬ 
ter a (sin 2 2'Olsnd = 4a 2 ). 

The most popular models which explain the solar, at¬ 
mospheric and LSND signatures (and the null results ob¬ 
tained by other experiments) via neutrino oscillations in¬ 
voke the existence of a sterile neutrino P s . One example 
of a (2 + 2) model is the following, which is taken from 

M: 


1 P s ^ 


O 

o 


( Vo \ 



-75 75 P P 


+1 



P 75 75 


Vi 

\ '+■) 


^ 0 0 -775 T75 ) 

\y 3 J 


( 12 ) 


where one pair of nearly degenerate mass eigenstates has 
maximal v e —> v s mixing for solar neutrinos and the other 
pair has maximal Vy —> v T mixing for atmospheric neutri¬ 
nos. Small inter-doublet mixings through the (3 parame¬ 
ter accomodates the LSND result (sin 2 'Olsnd = 8/3 2 ). 

Recent experimental results 0 show that pure v e —> 
v s solar oscillations are excluded at high significance. We 
therefore consider a more general (2 + 2) scenario, in 
which solar neutrinos can undergo any combination of 
v e —> v s and v e —> v T oscillations, while atmospheric 
neutrinos can undergo any combination of Vy —> v T and 
Vy —> v s oscillations. We follow the procedure in |l3] 
to obtain this more general mixing starting from Eq, 12, 
by substituting the ( P Sl P T ) states with the rotated states 
(*«) = 


The results for the Py,P T ,P e —> P e adiabatic oscillation 
probabilities, the permutation factor p , and the LSND 
oscillation amplitude sin 2 2$ as a function of the mixing 



cos <p s sin ip s 
- sin <p s cos p s 



(13) 
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Model 

Mixing 

P^— 

Pr—*e 

Pe—>e 

P 

sin 2 2i)lsnd 

Normal (1+1) 

Eq. 

io 


sin 2 ■& 

0 

cos 2 1 ? 

sin 2 & 

sin 2 2i? = 4p(l — p) 

LSND-inverted (1+1) 

Eq. 

10 


cos 2 $ 

0 

sin 2 1 ? 

cos 2 1 ? 

sin 2 2i? = 4p(l — p) 

Normal (2+1) 

Eq. 

11 


4 Lt 


1 

a 2 /(l + a 2 ) 

4 a 2 = 4p/(l — p) 

LSND-inverted (2+1) 

Eq. 

11 


1 

la 2 

12 

4 a 

(l + |a 2 )/(l + o 2 ) 

4a 2 = 4(1 — p)/(p — §) 

Normal (2+2) 

Eq. 

14 


P 2 

P 2 

1 

2 

4/17(1 + 4/3 2 ) 

8(3 2 = 2p/(l — p) 

LSND-inverted (2+2) 

Eq. 

14 


1 

2 

1 

2 

P 2 

1/(1+/3 2 ) 

1 

H 

OO 

II 

(N 

OO 

Normal (3+1) 

Eq. 

15 


7 2 

0 

1 

2 

2 7 2 /(l + 2 7 2 ) 

4 7 2 A 2 = 25 2 p/(l-p) 

LSND-inverted (3+1) 

Eq. 

15 


0 

1 

2 

7 2 

7(1 + 2 7 2 ) 

4 7 2 A 2 = 25 2 {l-p)/p 


TABLE II: Results on the probabilities for a to emerge from the SN as a i/ e , the permutation factor p of Eq.^, 

and the LSND oscillation amplitude sin 2 2i)lsnd, for the various neutrino mixing schemes considered. 


where the rotation angle ip s fixes the sterile component 
in the atmospheric doublet (0 < ip s < tt/2 ). Eq|l| then 
becomes: 


^ Vs ^ 


7 


7 


V 7 ) 



COS tp s 

COS (fis 

sin ip s 

sin ifs 

\ 

( 7 \ 

V2 

+2 

+2 

+2 



1 

n/2 

1 

P 

P 



7 

P 

-p 

1 

a/2 

1 

V2 



7 

sin 

sin <p s 

COS <£> s 

COS ip s 

) 


\ 7 / 

V2 

+2 

+2 

+2 



(14) 


which contains Eq.|l2| in the specific case ip s = 0. 
We note that the LSND oscillation amplitude formula 
sin 2 2 "Dlsnd = 8/? 2 holds also for the more general case 
of Eq.[l4], and that our results are independent of the 
value of tp s (see Tab. [n]). 

Another possible 4-neutrino model has a (3 + 1) hier¬ 
archy; as an example for this model, here we consider the 
following mixing, which is also taken from [p+|: 


^ 7 ^ 


(ii 2 A 0 

1 1 1 A 

7 


2 2 +2 ° 

7 


1 4 75 0 

V 7 ) 


U*"*7 1/ 


/ v x / 

7 

7 

V 7 / 

(15) 


where the solar and atmospheric oscillations are approxi¬ 
mately described by oscillations of three active neutrinos, 
and the LSND result by a coupling of 7 and v e through 
small mixings with D s that has a mass eigenvalue widely 
separated from the others (sin 2 2 Plsnd = 4^/ 2 <5 2 ). For 
the (3 + 1) scenario, the constraint given by the per¬ 
mutation probability p is not sufficient to determine the 
LSND oscillation amplitude sin 2 2 ’&lsnd- Therefore, the 
constraint on | L/^, 4 1 2 = S 2 given by the CDHS and Super- 
K experiments will also be used, as explained later. 

We should note that the mixing matrices defined in 
Eqs. [To|- [l5| are approximations in the sense that the ma¬ 
trices are unitary only up to order 0(a, (3, 7 , S). These 
are the parameters in the mixings responsible for LSND- 
type oscillations, which we let float for our analysis, but 
we know they are small. 


In order to determine the permutation factor for 
the mixing models, we also need to specify the neutrino 
mass hierarchy. In this paper, we consider for each mix¬ 
ing model both the cases of a “normal” and a “LSND- 
inverted” mass hierarchies. By “normal” hierarchy, here 
we mean that rrii > rrij for i > j, where to, is the 
mass eigenvalue for the |z/j) state. We define the “LSND- 
inverted” hierarchies as the ones obtained substituting 
Aitilsnd —> — A.rriLSND in the normal hierarchies, with¬ 
out changing the hierarchy of the eventual solar and at¬ 
mospheric splittings (see Tab. |lj); Aitilsnd is the neu¬ 
trino mass difference responsible for LSND oscillations. 

A common feature to all the mixing schemes is ap¬ 
parent in Tab. |n]. In the adiabatic approximation, nor¬ 
mal mass hierarchies predict small permutation factors, 
while an almost complete permutation would be present 
for LSND-inverted hierarchies. 

Given the specific neutrino mixing models considered 
here, we can now partially address the question whether 
the adiabatic approximation is applicable in this context. 
At a resonance, where the non-adiabaticity is maximal, 
this is a good approximation if the width of the resonance 
region is large compared with the local neutrino oscilla¬ 
tion length. The width of the resonance is, in turn, de¬ 
termined by the characteristic length scale of the radial 
matter density variations at the resonance. While there 
are reliable models for the matter density profile of the 
progenitor star, there are still uncertainties on the profile 
seen by neutrinos in their free-streaming propagation. 

It is now thought that neutrino heating of the proto¬ 
neutron star mantle drives the supernova explosion, 
which would happen with a ~ Is delay after the creation 
of the shock-wave, ultimately responsible for the explo¬ 
sion; during this delay, the shock-wave would be stalled at 
a radius of ~ 200 km from the neutron star, correspond¬ 
ing to a density p ~ 10 9 — 10 10 g/cm 3 [Q. Therefore, the 
density profile in the proximity of the stalled shock-wave, 
which is difficult to model reliably, is a potential site for 
non-adiabatic oscillations. 

In Fig. [i] we show the energy splittings between the lo¬ 
cal neutrino energy eigenvalues Ei, as a function of mat¬ 
ter density, for all eight neutrino models considered here. 
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°gio(p(g/cm 3 )) 


FIG. 1: Splittings between energy eigenvalues versus mat¬ 
ter density p for various neutrino mass and mixing mod¬ 
els. Solid, dashed, dotted lines show the splittings E 12 , E 23 , 
E 34 , respectively (see text). The local minima correspond to 
MSW-resonances. Model: a)normal (1+1); b)inverted (1+1); 
c)normal (2+1); d)inverted (2+1); e)normal (2+2); f)inverted 
(2+2); g)normal (3+1); h)inverted (3+1). Apart from the in¬ 
consequential <-> V r one in Figjljf, no MSW-resonances 
occur before the antineutrinos reach the stalled shock-wave 
(hatched area). 


For an n-neutrino model, we plot Ei^ + \ = £) — Ei + 1 , 
where i = 1 ,... ,n — 1; the eigenvalues are ordered such 
that Ei > E 2 >■■■> E n . Clearly, a resonance cor¬ 
responds to a local minimum in one of the curves. As 
can be seen from Fig. |l], all the resonances (except the 
inconsequential one in Fig.[l]f between and v T ©) 
lie at densities well below the stalled shock-wave den¬ 
sity of p ~ 10° — IO 10 g/cm 3 . Therefore, the impact of 
level crossing between propagation eigenstates is likely to 
be small even where the neutrinos encounter the shock- 
wave. 

If the SN neutrinos cross the Earth on their way to 
the detector, as for example happened for the SN1987A 
Ve's detected by the Kamiokande-2 and IMB-3 detectors, 
it is also necessary to evaluate the importance of Earth 
matter effects in the neutrino propagation. Clearly, for 


neutrino oscillation models where no solar splitting is in¬ 
volved (for example the (1 + 1) and (2 + 1) models in this 
paper), this effect is negligible. In the models where such 
a splitting is allowed ( i.e. the (2 + 2) and (3 + 1) mod¬ 
els considered here), the situation is more complicated. 
However, the Earth matter effects have been shown to be 
small in this case as well for a large fraction of the SN v e 
energy spectrum (below ~ 40 MeV) |[Uil . j2(j|| , and for the 
sake of simplicity will not be considered further. 


III. CONSTRAINTS ON LSND FROM SN1987A 
OBSERVATIONS 

Twenty v e events from the supernova SN1987A were 
observed by the Kamiokande-2 (Kam-2) and IMB-3 de¬ 
tectors. Kam-2 saw 12 events with an average energy 
of ( Edet) = 14.7 MeV, IMB-3 (which had a higher 
energy threshold than Kam-2) detected 8 events with 
{Edet) = 31.9 MeV 

From a comparison of the measured energy spectra 
(Fp e ) with theoretical models of neutrino emission (F 3 
and Fg ), it is possible to infer the permutation factor 
p in Eq. [I], SN1987A observations are consistent with 
no-oscillations {i.e. p = 0). In Appendix A, we derive a 
conservative upper bound onp oip < 0.22 at 99% CL, by 
applying a Kholmogorov-Smirnov test on the joint Kam- 
IMB dataset and a range of supernova neutrino emission 
models. 

One important result of our analysis is immediately 
apparent from the values of the permutation factors p 
as a function of the mixing parameters in Tab. [U|. and 
from the fact that the value of p inferred from SN1987A 
data has to be less than 0.22 at 99% CL. The four mix¬ 
ing schemes considered, explaining the LSND effect via 
a LSND-inverted neutrino mass hierarchy, are all incom¬ 
patible with SN1987A data. 

We now consider the normal hierarchy cases. For the 
(1 + TL (2 + 1) and (2 + 2) models with the mixings of 
Eqs. 1(1112, the bound on the permutation factor p un- 
ambigously determines the constraint on the LSND os¬ 
cillation amplitude sin 2 2 ’Olsnd (see Tab. At 99% 
CL, SN1987A data provide no constraints on the (2 + 1) 
model, and a constraint which is weaker than existing 
bounds from the accelerator ex per iment Karmen p3[ and 
the reactor experiment Bugey [2^] for the (1 +1) and 
(2 + 2) models (see Fig. [|r). Therefore, these models are 
compatible with the SN1987A data. 

As already mentioned, for the (3 + 1) model, the per¬ 
mutation factor does not fully determine the LSND os¬ 
cillation amplitude: sin 2 2 'Olsnd depends not only on 
p, but also on | C /^4 1 2 = S 2 . Here we use the Am 2 LSND - 
dependent constraints on S 2 from the ^-disappearance 
experiments CDHS (for Am 2 Q Mn > 0.3 eV 2 ) and Super- 
K (for Am 2 LSND < 0.3 eV 2 ) J25|. Moreover, another 
complication arises in evaluating exclusion regions for 
(3+1) models: given the 99% CL upper bounds on 
7 2 = |C/ e 4 1 2 from SN1987A and S 2 = \U 2 A \ from CDHS 
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FIG. 2: 99% CL LSND allowed region [jj and 99% CL exclusion regions for the neutrino mixing schemes considered in the 
text and with normal mass hierarchy. The exclusion regions are estimated as in a) shows the exclusion regions for the 

(1 + 1), (2+1) and (2 + 2) models, b) for the (3 + 1) model. The exclusion regions refer to experimental data from the following 
experiments, a) Dotted line: Karmen; dashed line: Bugey; dark solid line: SN1987A for the (2 + 2) model; light solid line: 
SN1987 for the (1 + 1) model; SN1987A data provide no constraints at 99% CL for the (2 + 1) model, b) Dotted line: Karmen; 
dashed line: Bugey, CDHS and Super-K; solid line: SN1987A, CDHS and Super-K. 


Model 

SN1987A constraint on 
LSND region (99% CL) 

Normal (1 + 1) 

partially excluded (Fig. ^,) 

LSND-inverted (1 + 1) 

excluded 

Normal (2 + 1) 

unconstrained 

LSND-inverted (2 + 1) 

excluded 

Normal (2 + 2) 

partially excluded (Fig. ^|a) 

LSND-inverted (2 + 2) 

excluded 

Normal (3 + 1) 

partially excluded (Fig. ^b) 

LSND-inverted (3 + 1) 

excluded 


TABLE III: Summary of the SN1987A constraints on the 
LSND allowed region, for the various models considered in 
this paper; see Fig. |i| also. 


and Super-K, what is the 99% CL upper bound on 
sin 2 2'Olsnd = 4+ 2 <5 2 ? We follow the method described 
in @ to estimate this bound. The same method is ap¬ 
plied to estimate the 99% CL upper limit on sin 2 2 Alsnd 
coming from Bugey (for + 2 ) and CDHS and Super-K (for 
<5 2 ), that is without using the SN1987A data. The re¬ 
sults for the (3+1) model with normal neutrino mass 
hierarchy and mixing given by Eq. |l5]is shown in Fig. j|b. 
Also for this model, we find that existing constraints (the 
Bugey constraint on S 2 , in this case) are stronger than 
the SN1987A one. 

Tab. Ill summarizes the SN1987A constraints obtained 
in this paper on the LSND allowed region, for the various 
neutrino mass and mixing models considered. 


IV. CONCLUSIONS 


We have investigated the effect that 3- and 4-neutrino 
oscillation schemes would have in modifying the energy 
spectrum of supernova P e ’s. Throughout the paper, we 
apply the adiabatic approximation for the antineutrino 
propagation in the supernova environment and neglect 
Earth matter effects. Moreover, we have used our results 
to test the compatibility between the SN1987A data and 
the LSND evidence for 7+ —» v e oscillations. 

We have provided specific relations for the permutation 
factor, which gives the admixture of a higher energy flux 
to the original v e flux at production from v^.,u T —» v & os¬ 
cillations, for various neutrino mass and mixing models. 
The permutation factor may be measurable with good 
accuracy in future supernova experiments. 

Based on SN1987A data only, which seem to indicate 
a small (if non-zero) value for the permutation factor, 
we are able to exclude all of the four models consid¬ 
ered which would explain the LSND signal via a “LSND- 
inverted” neutrino mass hierarchy, as defined in the 
text. For the normal mass hierarchy schemes considered, 
SN1987A data do not provide any stronger constraints 
on the LSND allowed region for oscillations than those 
already obtained with reactor, accelerator and atmo¬ 
spheric neutrinos; additional experimental input is nec¬ 
essary to unambiguously discern the neutrino mass and 
mixing properties. Undoubtedly, the detection of super¬ 
nova neutrinos by present or near-term experiments |27jj 
would prove very useful in this respect. 
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APPENDIX A: UPPER BOUNDS ON THE 
PERMUTATION FACTOR FROM SN1987A DATA 

In this appendix, we discuss the statistical methodol¬ 
ogy and the physics assumptions used to estimate the 


upper bound on the permutation factor p quoted in the 
text, p < 0.22 at 99% CL. We use the same statistical 
methodology as in |2^], that is we use the Kholmogorov- 
Smirnov test on the joint Kam-IMB dataset to derive the 
upper bound. Most of the physics assumptions are iden¬ 
tical to those in p8[ . 

The expected energy spectrum for the positrons, ob¬ 
served in the Kamiokande and IMB detectors via the 
reaction v e p —> e + n, is: 


Bi(Edet) 


n pA r 
4 nD 2 J 0 


dE + Pi(Edet, E + )po } i(E + )av e p(E + + Q)F„ e (E + + Q) 


(Al) 


where i refers to either Kam or IMB, N Pi i is the number 
of target protons in the detectors, D the distance between 
the Large Magellanic Cloud and the Earth, E d et ( E + ) is 
the detected (true) positron energy, Q = m n — m p = 
1.29MeV ~ E Pe - E + , Pi(E det ,E + ) and rj 0 ,i(E + ) the 
energy resolution functions and efficiency curves taken 
from ^8|, ap eP (E + + Q) oc E+ the neutrino interaction 
cross-section taken from || (neglecting nuclear recoil), 
and finally Fp e (E + +Q ) the neutrino flux at the detector 
taken from Eqffi. We assume “unpinched” Fermi-Dirac 
distributions for the fluxes F p , a = e, /x, appearing in 
Eq.|: 


FUE)o c 


E 2 


(F p jT3(e E /^ + 1) 


(A2) 


where {E Pa ) — 3.15T a at the denominator ensures energy 
equipartition. 

The cumulative distribution function used for the 
Kholmogorov-Smirnov test is: 

E(Edet) = / dE(n Ka m(E)+n IM B(E)) (A3) 

Jo 

Fig.|] shows the upper bound on the permutation fac¬ 
tor p obtained from SN1987A data, at 99% CL, as a 
function of the average energies (E Pe ), (E„ ). As ex¬ 
pected, the bound becomes more stringent for supernova 
models in which the neutrino average energies are higher. 
SN1987A data are incompatible at 99% CL with all su¬ 
pernova neutrino models predicting (E Pe ) > 16.6 MeV, 
for all values of p and {E Dfj ). We adopt a conserva¬ 
tive approach, and quote as the upper bound on p the 
largest value for supernova neutrino models in the range 


14 < (E Pe ) < 17 MeV, 23 < {E P J < 27 MeV, that is the 
one corresponding to (E Pe ) = 14 MeV, (E^) = 23 MeV 
(cross in Fig-D- 



FIG. 3: Solid lines: isocontours for the upper bounds on 
the permutation factor p at 99% CL obtained from SN1987A 
data, as a function of the v e and v p average energies predicted 
at production by supernova models; rectangle with dashed 
border: range of energies allowed by present models; cross: 
model chosen to derive the (conservative) upper bound on p 
used in the text. The region (E Pe ) > 16.6 MeV is excluded 
at 99% CL for all values of p. 
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